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Abstract.  A  reformulation  of  the  planetary  geostrophic  equations  (PGEs)  with  inviscid  balance  equation 
is  proposed  and  the  existence  of  global  weak  solutions  is  established,  provided  that  the  mechanical  forcing 
satisfies  an  integral  constraint.  There  is  only  one  prognostic  equation  for  the  temperature  field  and  the  ve¬ 
locity  field  is  statically  determined  by  the  planetary  geostrophic  balance  combined  with  the  incompressibility 
condition.  Furthermore,  the  velocity  profile  can  be  accurately  represented  as  a  functional  of  the  temperature 
gradient.  In  particular,  the  vertical  velocity  depends  only  on  the  first  order  derivative  of  the  temperature. 
As  a  result,  the  bound  for  the  L°°(0,  L2)  DL2(0,  t\\  H1)  norm  of  the  temperature  field  is  sufficient  to  show 

the  existence  of  the  weak  solution. 
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1.  Introduction 

The  planetary  geostrophic  equations  (PGEs)  have  played  an  important  role  in  large-scale  ocean  circulation 
since  the  pioneering  work  of  A.  Robinson  and  H.  Stommel  [6]  and  P.  Welander  [12].  This  system  arises  as 
an  asymptotic  approximation  to  the  primitive  equations  (PEs)  for  planetary-scale  motions  in  the  limit  of 
small  Rossby  number.  The  PGEs  are  considerably  simpler  than  the  PEs,  but  retain  the  dynamics  necessary 
to  represent  the  large-scale,  low-frequency  dynamics  of  the  mid-latitude  oceans. 

The  PGEs  with  viscous  geostrophic  balance  has  been  analyzed  at  the  PDE  level  in  recent  articles.  See 
[1,  8,  9]  for  relevant  discussions.  One  distinguishing  feature  of  the  PGEs  is  that  there  is  only  one  prognostic 
equation  in  the  system  for  the  temperature  field;  the  velocity  field  is  diagnostically  determined  by  the 
planetary  geostrophic  balance.  The  addition  of  a  diffusion  term  in  the  geostrophic  balance  equation  is  for 
the  sake  of  simplicity  in  mathematical  analysis,  due  to  the  lack  of  regularity  for  the  velocity  field  by  a 
straightforward  manipulation . 

In  this  article,  we  consider  the  original  formulation  of  the  PGEs,  with  no  viscous  term  in  the  geostrophic 
balance  equations.  In  such  a  formulation,  both  horizontal  and  vertical  velocity  profiles  are  accurately  rep¬ 
resented  as  functionals  of  the  temperature  gradient.  The  representation  formula  for  the  horizontal  velocity 
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field  is  based  on  the  planetary  geostrophic  balance.  Since  every  variable  can  be  uniquely  determined  by  the 
combination  of  its  mean  (in  vertical  direction)  and  its  vertical  derivative,  the  horizontal  velocity  turns  out  to 
be  the  solution  of  a  differential  equation  at  each  fixed  horizontal  point,  depending  only  on  the  temperature 
gradient.  The  vertical  velocity  can  be  recovered  by  the  continuity  equation.  Using  the  special  form  of  the 
Coriolis  parameter,  we  arrive  at  a  two-point  boundary  value  ordinary  differential  equations  (O.D.E.)  in  the 
vertical  direction  at  each  fixed  horizontal  point  for  the  vertical  velocity,  with  the  right  hand  side  depending 
only  on  the  first  order  derivative  of  the  temperature  field. 

The  new  formulation  is  derived  in  Section  2  and  the  existence  of  the  global  (in  time)  weak  solution  of 
the  reformulated  PGEs  is  provided  in  Section  3.  The  approach  of  Galerkin  approximation  is  used.  Standard 
energy  estimate  for  the  temperature  equation  gives  the  bound  of  the  -L°°(0,  A;  L2)  fl  L2(0,  A;  R1)  norm  of 
the  temperature  variable,  which  in  turn  shows  the  bound  of  the  L2( 0,  t\  \  L2)  norm  of  the  horizontal  velocity. 
Moreover,  the  L2(0,  t\\  L2)  norm  of  the  vertical  velocity  is  also  uniformly  bounded  since  it  satisfies  the  second 
order  O.D.E.,  in  which  the  force  term  only  involves  the  temperature  gradient.  The  compactness  for  the  time 
derivative  of  the  temperature  field  can  be  established  in  a  similar  manner.  Thus  the  existence  of  the  global 
weak  solution  is  proven. 

2.  Reformulation  of  the  Inviscid  Planetary  Geostrophic  Equations 

The  non-dimensional  PGEs  can  be  written  as 

r,  +  („.v)r  +  «f  =  (^a  +  T_82)t, 

fk  x  v  +  Vp  =  F , 
dp  T 

dz 

V-t?  +  dzw  =  0  , 

where  T  represents  the  temperature,  v  =  (u,  v)  the  horizontal  velocity,  w  the  vertical  velocity,  and  p  the 
pressure.  The  term  fkxv  corresponds  to  the  Coriolis  force  with  /  depending  only  on  the  latitude  y.  As 
a  typical  example  used  in  geophysical  literatures,  its  f3— plane  approximation  is  given  by  f  =  fo  +  (3y ■  The 
parameters  1/Rti,  1/Rt2  stand  for  the  horizontal  and  vertical  heat  conductivity  coefficients.  The  operators 
V,  V^,  V-,  A  stand  for  the  gradient,  perpendicular  gradient,  divergence  and  Laplacian  in  horizontal  plane, 
respectively.  For  simplicity,  we  set  Ki  =  1/Rt-i,  n2  —  l/Rt2-  The  forcing  term  F  =  ( Fx,Fy )-L  appearing  in 
the  geostrophic  balance  equation  (2.1)2  comes  from  the  wind  stress  at  the  ocean  surface,  which  is  a  boundary 


layer  approximation.  It  may  or  may  not  depend  on  the  vertical  variable  z.  For  simplicity,  we  assume  in  this 
article  F  =  F(x,y)  =  (^Fx(x,y),Fv(x,y)^J.  The  discussion  of  a  general  case  can  be  carried  out  in  the  same 
fashion  and  does  not  add  any  mathematical  difficulty.  See  the  relevant  references  on  both  the  physical  and 
mathematical  descriptions  of  the  PGEs  in  [2,  3,  4,  5,  6,  8,  9,  10,  12],  etc. 

The  computational  domain  is  taken  as  M  =  A4o  x  [— Hq,0],  where  A4o  is  the  surface  of  the  ocean.  The 
boundary  condition  at  the  top  and  bottom  surfaces  are  given  by 


(2.2) 


w  =  0 
w  =  0 


and  H2  , 

az 

,  dT  n 

az 


at  z  =  0 , 
at  z=  —H0 , 
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where  the  term  T*  represents  the  heat  flux  at  the  surface  of  the  ocean.  Usually  T*  can  be  taken  as  either 
a  fixed  heat  flux  function  or  of  the  form  T *  =  —  a(T  —  9*),  where  9*  is  a  reference  temperature.  Both 
boundary  conditions  can  be  dealt  with  in  an  efficient  way.  In  this  article  for  simplicity  we  choose  T f  as  a 
given  flux.  On  the  lateral  boundary  section  dM.$  x  [— Hq,  0],  the  fixed  boundary  condition  is  prescribed  for 
the  temperature  field 

(2.3)  T  =  Tlb,  onSM0x[-ff0,0], 

where  Tib  is  a  given  distribution.  As  will  be  shown  later,  the  purpose  of  the  choice  of  a  Dirichlet  boundary 
condition  is  to  facilitate  the  analysis  of  the  system  at  the  PDE  level,  although  the  no-flux  boundary  condition 
for  the  temperature  field  is  physically  more  relevant.  The  boundary  condition  (2.3)  can  also  be  viewed  as  an 
approximation  such  that  the  disturbance  of  oceanic  circulation  motion  is  far  away  from  the  lateral  boundary. 
For  simplicity  of  the  presentation,  we  set  the  homogeneous  profile  Tib  =  0  in  the  theoretical  and  numerical 
analysis.  There  is  no  real  change  for  the  non-homogeneous  case.  The  normal  component  of  the  vertically 
averaged  horizontal  velocity  turns  out  to  have  a  vanishing  flux 

(2.4)  v  n  =  0,  on  dMo , 


which  is  compatible  with  the  continuity  equation  (2.1)4.  We  recall  that  the  average  (in  the  vertical  direction) 
of  any  3-D  field  g  is  given  by  g(x,  y)  =  S-h0  2/>  z)  ^ 2 •  See  [4,  8]  for  a  detailed  explanation  for  the 

choice  of  this  nonlocal  boundary  condition  in  the  case  where  no  viscosity  is  present  in  the  geostrophic  balance 
equation. 

We  now  derive  an  equivalent  formulation  of  the  system  of  PGEs  (2.1)-(2.4).  The  key  point  in  this 
reformulation  is  that  both  the  horizontal  and  vertical  velocity  variables  can  be  determined  by  the  first  order 
derivative  of  the  temperature  field.  This  makes  valid  the  analysis  of  the  well-posedness  for  the  system. 

The  horizontal  velocity  field  is  the  solution  of  the  following  system 


(2.5) 


*  -Tv 
dzu  =  —jL 


T 

°-v  =  7’ 


,  _  -dyF*  +  dxFy 

u(x,y)=ue,  v{x,y)  =  ve  =  - 


9yf 


where  ue  is  explicitly  given  below. 

Equation  (2.5)  is  derived  from  the  geostrophic  equation  and  hydrostatic  equation.  Taking  the  vertical 
derivative  of  the  geostrophic  balance  equation  fkxv  +  Vp  =  F  gives  the  thermal  wind  equation 


(2.6) 


v,  = 


V±Pz  +  dzFA 


f 


\7±T 

7 


i.e. 


Zlv 

f 


If 

f  ’ 


where  the  hydrostatic  balance  dp/dz  =  T  and  the  independence  on  z  of  F  and  /  =  f(y)  were  used.  In  other 
words,  the  profile  vz  can  be  expressed  by  the  temperature  gradient. 

Meanwhile,  integrating  the  geostrophic  balance  equation  fkxv  +  Vp  =  F  in  the  vertical  direction  and 
dividing  by  FIq ,  we  find 


(2.7) 


fk  x  v  +  Vp  =  F . 
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Applying  the  curl  operator  V-1  to  (2.7)  results  in 

(2.8)  (dvf)v  +  f(vy  +  ux)  =  V  x  F  =  -dvFx  +  dxF»  . 

Moreover,  the  continuity  equation  V  •  v  +  dzw  =  0  and  the  boundary  condition  for  the  vertical  velocity 
w(  •,  0)  =  tu(  ■,  —Ho)  =  0  show  that  the  averaged  horizontal  velocity  field  is  divergence- free,  i.e., 

(2.9)  V  •  tJ  =  0  . 


The  combination  of  (2.8)  and  (2.9)  yields 

(2.10) 


_  ,  -dyFx  +  dxFy 

ve  =  v{x,y)  =  - 


dyf 


By  taking  the  tangential  part  of  (2.7)  and  applying  (2.4),  one  obtains  the  vertically  averaged  tangential 
pressure  gradient  (and  thus  the  pressure,  aside  from  an  arbitrary  constant)  on  the  boundary  from  the 
tangential  component  of  the  forcing  F: 


(2.11) 


—  =  F  ■  r 
Or  ’ 


where  r  is  the  unit  tangential  vector  on  the  boundary.  Integrating  this  relation  around  the  boundary,  one 
obtains  the  constraint  that  the  line  integral  of  F  around  the  boundary  must  be  zero: 


(2.12) 


F  ■  t  dl  =  0  . 


IdM, 


This  means  that  the  forcing  must  not  give  a  net  torque  on  the  fluid. 

From  (2.9),  we  can  find  a  2-D  mean  stream  function  ip(x ,  y)  for  the  vertically  averaged  velocity  field,  such 
that  (u.  v)  =  (—dyip,dxip).  Moreover,  the  boundary  condition  (2.4)  indicates  that  ip  is  a,  constant  on  the 
lateral  boundary.  For  simplicity  of  the  discussion,  we  take  ip  =  0  on  <9A4o-  In  addition,  we  denote  71(2/0), 
72(2/0)  by  the  ^-coordinates  of  the  intersection  points  between  <9A4o  and  y  =  j/o-  The  mean  stream  function 
ip  and  the  mean  velocity  u  can  be  determined  by  the  kinematic  relationship  and  formula  (2.10): 

VxF 


(2.13) 


1pe(x>y)  = 


71  (y)  dyf 


■  dx' ,  ue(x,y)  =  —dyipe(x,y) , 


with  71(2/)  being  a  point  on  <9A4o-  Evaluating  ip  at  another  boundary  point  (72(2/),  y)  with  the  same  y- value, 
we  obtain  an  additional  constraint  on  the  forcing: 

C72(y)  y  x  F 

J’niv) 


(2.14) 


dvf 


dx'  =  0 , 


since  ip  is  identically  0  on  the  lateral  boundary.  Constraint  (2.14)  amounts  to  saying  that  the  average  forcing 
across  the  domain  at  a  fixed  y  must  not  give  a  torque  on  the  fluid. 

The  combination  of  (2.6),  (2.10)  and  (2.13)  leads  to  the  system  (2.5).  By  the  representation  formula 
valid  for  any  3-D  variable  g: 

(2.15)  g(x,y,z)=  f  gz(x,y,z1)dz1  +  g(x,y)-^~  [  [  gz{x,y,  zi)  dz-i  dz , 

J-H0  U0  J  —  Ho  J~H0 
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the  solution  of  (2.5)  can  be  expressed  explicitly  using  an  integration  formula: 

rz  t  1  fz  T. 

(2.16a)  u(x,  y,z)  =  —  /  (x,y,  z,)dz,  +ue(x,y)  +  —  /  /  -f(x,  y,  z,)  dz,  dz  , 

J  —  Hq  J  n0  J  —  H0  J-Hq  I 

fz  Tx  1  [z  Tx 

(2.16b)  v(x,y,z)  =  (x,  y,  z,)  dz,  +  ve(x,  y)  -  —  /  /  ~jr(x,y,  zx)  dzi  dz , 

J-H0  J  UQ  J-H0  J-H0  J 

with  ue,  ve  given  by  (2.10)  and  (2.13). 

The  vertical  velocity  can  be  calculated  by  integrating  the  horizontal  divergence  of  the  horizontal  velocity 
field,  due  to  the  incompressibility  condition 


(2.17) 


w{x,  y,z)  =  -  V-v(a %y,  s)  ds . 

J-H0 


The  substitution  of  (2.16)  into  (2.17)  gives 


(2.18)  w(x,y,z)=  [  [2  ^j^(x,y,z1)dz1dz2-^-(z  +  H0)  [  [  {-^j^(x,y,  Zl)  dz,  dz . 

J  —  Ho  J  —  Ho  J  U0  J  —  H0  J-H0  J 

Therefore,  w  can  be  expressed  as  a  functional  of  the  temperature  gradient,  like  the  horizontal  velocity  v. 

The  vertical  velocity  can  also  be  represented  as  the  solution  of  a  differential  equation.  By  taking  the 
vertical  derivative  of  the  continuity  equation 

V  •  vz  +  <9? u;  =  0  , 


(2.19) 

combined  with  (2.6),  we  arrive  at 

(2.20)  d2zw  =  -dx(uz)  -  dy{vz )  =  -  ^(y) 


Tx\  _  ( dyf)Tx 

f>  f2 


It  can  be  observed  that  the  second  order  derivatives  for  the  temperature  field  cancel  each  other  due  to  the 
special  form  of  the  Coriolis  parameter  /  =  f(y).  Therefore,  the  vertical  velocity  w  can  be  reformulated  as 
the  solution  of  the  following  system  of  second  order  O.D.E.s 


,2„„  _  (dyf)Tx 


dzw  = 


(2.21) 


/ 2 


w  =  0  ,  at  2  =  0,  —  Hq  , 


in  which  the  right  hand  side  includes  only  the  first  order  derivative  of  the  temperature.  This  key  point  is 
crucial  to  the  analysis  presented  below. 

We  then  have  the  following  formulation,  where  the  velocities  are  expressed  as  functionals  of  the  temper¬ 
ature  gradient. 

Temperature  Transport  Equation 


(2.22a) 


dT  /  \ 

Tt  +  (vX7)T  +  w—  =  (yA  +  K23?jT, 


dT 

dz 


=  Tf  ,  at  2  =  0, 


dT 

dz 


=  0 ,  at  2  =  —  H0 


T  =  0,  on  dM0  x  [—#0,0]; 
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Recovery  of  the  horizontal  velocity 


(2.22b) 


d~u  = 


Ty 

f  ’ 


u(x,  y)  =ue, 


v(x,y)  =  ve ; 


Recovery  of  the  vertical  velocity 


(2.22c) 


d'i'W  = 


0 dyf)Tx 

f2 


w  =  0  ,  at  z  =  0  ,  —  H0  ■ 


Remark  2.1.  It  is  observed  that  the  alternate  formulation  (2.22a-c)  is  equivalent  to  the  original  formulation 
(2.1)-(2.f)  of  the  PGEs,  from  which  they  were  derived.  Indeed,  to  recover  (2.1)-(2.f)  from  (2.22a-c),  we 
need  to  show  that  f>  =  V  •  v  +  dzw  =  0.  A  simple  calculation  utilizing  (2.22b)  and  (2.22c)  leads  to 

(2.23)  dz(j)  =  dz(\7  ■  v  +  dzw^j  =  dx(dzu )  +  dy(dzv )  +  d2w  =  -dx(y-)  +  dv(~j)  + 


(2.24)  (f>  =  •  v  +  dziv'j  =  V  •  v  +  dzw  =  0  , 

since  the  average  of  v  is  divergence- free  in  the  horizontal  plane.  The  combination  of  (2.23)  and  (2.2 f)  results 
in  the  incompressibility  condition.  In  addition,  a  direct  calculation 

(2.25)  V  x(fkxv-F)  =  (dvf)v-V  xF  =  0, 


indicates  the  existence  of  a  mean  pressure  field  p  such  that  fkxv  +  Vp  =  F.  Accordingly,  we  define  the 
total  pressure  field  as 


(2.26) 


P(x,y,z)=  [  T(x,y,s)ds+p(x,y)  -  f  f  T(x,  y,s)  dsdz . 
J-Hn  uO  J-Hn  J-Hn 


Clearly,  the  hydrostatic  balance  is  satisfied  by  taking  the  vertical  derivative  of  (2.26).  The  geostrophic  balance 
equation  can  also  be  verified  by  using  the  integration  formulas  for  v  in  (2.16),  which  comes  from  the  recovery 
equation  (2.22b),  combined  with  the  horizontal  gradient  of  (2.26): 

(2.27)  Vp=  [  VTds  +  Vp--^  [  [  S/Tdsdz. 

J  —  H0  ho  J-h0  J-h0 


All  the  boundary  conditions  presented  in  (2.1)-(2.f)  are  included  in  the  system  (2.22).  Hence,  ( T,u,p )  is 
also  a  solution  of  (2.1)-(2.f).  This  completes  the  proof  of  the  formal  equivalence  of  smooth  solutions  between 
the  two  formulations. 
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(3.1) 


3.  Existence  of  a  Global  Weak  Solution 

Before  starting  the  discussion  on  the  weak  solutions,  we  introduce  the  following  functional  setting: 
H  =  L2(M) ,  V  =  the  closure  of  Cf£t  0(M)  in  H1(M)  , 

=  {T  G  Cm(M)  |  T  =  0  ondMo*  [~H0, 0]}  . 


Note  that  the  introduction  of  C'[)lt  0(At)  is  motivated  by  the  boundary  condition  for  T  on  the  lateral  boundary 
sections.  Let  ( •,  •)  be  the  inner  product  in  L2(M)  =  L2(Mo  x  [—Ho,  0]),  and  ||-||  the  corresponding  L 2  norm. 
For  any  positive  final  time  t\  >  0,  the  functions  (T,  u,  v,  w) 

T  G  L°°(o,ti  ;#)  nL2(0,fi;E)  ,  u,v,w  G  L2  (o,h;  L2(M)^j  , 
are  called  a  weak  solution  of  the  original  PGEs  formulated  in  (2.22)  if 

| ' dt(T<j) )  +  uT(f>x  +  vT<j)y  +  wT(j>z  +  +  K2(dzT)  ■  (dz<j>)  ^  dxdz 

+n2  [  -(T(x,O)-0*)<f>(x,O)dx  =  O,  \/(j>  G  Cjat  0(A1)  D  H3(A4) , 

JMn  K2 


(3.2a) 


)M o  J-H0 


where 


(3.2b) 


rz  t  i  [0  rz  t 

i(x,y,z)  =  -  -jr(x,y,  z^dzi  +  ue(x,y)  +  —  /  /  -^-(x,y,z1)dz1dz 

J  —  Ho  J  J  —  Hq  J-H0  J 

>{x,y,z)=  [  -jr{x,y,z1)dz1+ve{x,y)  - [  [  -^-(x,y,z!)dzidz, 

J-H0  J  J-H0  J-H0  J 


(3.2c)  w(x,y,z)  =  fZ  f2  {d^(x,y,z1)dz1dz2-^(z  +  H0)  [  f  {M^{x,y,z1)  dzx  dz . 
J-H0  J ~Hq  J  J-Ho  J-Ho  J 

Theorem  3.1  Suppose  F  G  H2(Mq)  is  given  and  the  constraint  (2.14)  is  satisfied  so  that  ue,  ve  can  be 
consistently  determined.  Let  Tq  =  T(-,  0)  G  L2(M).  Then  there  exists  at  least  one  global  weak  solution  for 
the  PGEs  (2.22),  such  that  for  any  t\  >  0 


(3.3) 


(3.4) 


Tg  ni2(o,fi;E)  ,  dtT  G  L*  (o,  h  ;  H~2{M))  , 

u,v  ,w  G  L2(o,ti',  L2(M ) j  ,  dzu , dzv ,  d2w  G  L2  ^0,  ti ;  L2(M )^  , 

||T(-,f)||2  +  2^  (ki||VT(-,  s) || 2  +  K2||3ZT(-,  s)||2^  ds  <  1 1 To 1 1 2  +  C*t ,  for  0  <  t  <  h  , 

with  C*  =  a  f  (6*)2  dx  , 

Jm„ 


(3.5) 


I -Ho 


v  dz  =  0  ,  in  the  sense  of  distribution  , 
(0 

/  vdz-n  =  0,  on  dMo , 

I -h0 

w  =  0 ,  at  x  =  0 ,  —Hq  . 
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Proof.  The  proof  can  be  accomplished  by  the  Galerkin  procedure.  A  standard  energy  estimate  is  used  to 
obtain  the  uniform  bound  for  the  L°°(0,  fi ;  L2{ M.))  and  L2(0,ti ;  Hl{JV[))  norms  of  the  temperature  field. 


Let  {$j}  ->:l  C  H2(M)  be  the  eigenvectors  of  the  diffusion  operator  corresponding  to  the  eigenvalues  Ay, 
j  =  1,  2, ...,  such  that 


f  A  +  k282^J  <f >j  =  X j$j ,  Xj  — >  oo  , 


(3.6) 


<)'\>,  <l> , 

— —  =  —a—  ,  at  z  =  0  , 

OZ  H2 

=  0 ,  on  8M  o  x  [—Hq,  0] . 


=  0 ,  at  z  =  -H0  , 

OZ 


The  diffusion  operator  A  =  n±A  +  n2d2  with  the  given  boundary  condition  is  a  self-adjoint  linear  operator 
and  admits  a  compact  inverse.  Then  {<&„}n>1  defines  a  complete  orthogonal  basis  in  L2(M).  To  seek  a 
weak  solution  of  the  reformulated  PGEs  defined  in  (3.2),  we  find  an  approximate  solution  {Tm}  such  that 


(3.7) 


Tm{x,z ;  t)  =  '^2(3Jl(t)<$>j(x,z) 
i= i 


(3.8) 


dt 


(Pmi'&j'j  +  M  ^  VTm ,  V  $  j  ^  -\-  K2^dzTm,dz<& 

+k2  [  -(Tm(x,0)-0*)$j(x,0)dx  =  0,  j  =  1, 2, ...,  m , 

J  Mn  K2 


(3.9) 


( x,y,z)  =  -f  dAJ^{x,y,z1)dz1  +  ue{x,y)  +  f  f  ^^A(x,y,z1)dz1dz, 
J-Hn  J  n0  J-Hn  J-Hn  J 


Jy-L  m  i 

I -H0  J -H0  f 

i{x,y,z)=  [  ^AA(x,y,z-s)dzi  +  ve(x,y)  -  (x,  y,  zx)  dzi  dz , 

J-Hn  J  n0  J  —  Ho  J  —  H0 


I -H0  f 
8  T 

y~,x-Lm 


Jx±  m  , 

T 


(3.10) 

wm(x,y,z)  =  / 

J-H0  J-Ho 


22  (, dyf)dXTn 


(  V  ,  ,  1  ,  ^  „  s  r  fz  ( dyf)dxTm 

[x,y,z1)dz1dz2- —{z  + H0)  /  /  - (x,  y,  zx)  dzx  dz  , 

Ho  J-Ho  J-Ho  J 


(3.11)  Tm  |t=0 —  P-rnTo  , 

where  Pm  is  the  orthogonal  projection  operator  in  L2(M):  Pm:  L2(M )  — >  Span{$i,  ...,<f>m}. 

The  scheme  (3.8)  and  (3.11)  proposes  an  initial  value  problem  for  a  system  of  m  O.D.E.s,  with  the 
velocities  determined  by  (3.9)  and  (3.10).  Therefore,  it  is  straightforward  to  conclude  the  local  (in  time) 
existence  of  the  approximate  solution.  To  get  the  global  (in  time)  solution,  the  energy  estimates  are  necessary. 
We  observe  that  the  approximated  velocity  field  um  =  (vm,wm)  satisfies 

^  '  Vjn  T  OzWm  —  0  , 

(3.12) 

Vm  *  XL  l^^vio  X  [— //o,0]  Ve  '  XL  IcJATq  X  [— //o,0]  —  0  5  LUm  \z=0,  —  Hq  0  , 
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which  comes  from  the  construction  (3.9),  (3.10)  and  the  homogeneous  Dirichlet  boundary  condition  for  Tn 
on  the  lateral  boundary.  As  a  result,  we  find  by  integration  by  parts  that 

(3.13)  (vmTm,  VTm)  +  (wmTm,  dzTm^j  =  0  . 

Multiplying  (3.8)  by  and  adding  up  the  resulting  equations  leads  to 

~JTmf  +  Ki\\VTm\\2  +  K2\\dzTm\\2  =  -K2  \  —  (Tm(-,0)-e*)Tm(;0)dx 

2  at  JMo  a2\  / 

=  -a  f  (rm(;0 )-e*^Tm(-,0)dx 


a 

< - 

2 


a 

<  - 
~  2 


}  [  T2(,o)^  +  f  f  (e*fdx 

-  JM0  z  JMo 

f  (9*)2  dx . 

J  Mn 


Applying  the  Gronwall  inequality  to  (3.14)  results  in 


(3.15)  ||Tm(-,f)||2  +  2  [  MVTJ|2H-K2||dzTm||Vs<  ||r0||2  +  C*t,  with  C*  =  a  f  (9*)2  dx , 

Jo  J  Mn 


which  in  turn  indicates  that 


Tm  G  a  bounded  set  of  L°°  ^0,fi ;  L2(M)^j  D  L2  ^0,  ti ;  . 


Moreover,  by  the  recovery  formulation  (3.9),  (3.10),  we  have 

um  ,  i>m  ,  wm  G  a  bounded  set  of  L 2  ( 0,  fi ;  L2(M)  )  , 

(3.17)  .  K  (  ■  \ 

dzum  ,  dzvm  ,  d2wm  G  a  bounded  set  of  L2  (  0,  t\ ;  L2(M) )  . 

Furthermore,  we  need  an  estimate  for  <9tTm  so  that  we  can  obtain  compactness  and  a  strong  convergence 
result.  Consider  T  G  H2(M)  given  by  T  =  HjLi  Pj^j-  Equation  (3.8)  shows  that 

(c>tTm,f)=  ( 'dtTm,PmT ) 

13. 181  =  (vmTm^PmT)  +  (wmTm,dzPmf^  +  +  K2(dzTrn,dzPmf] 


JMo  K2 

Regarding  the  nonlinear  term,  we  have 


{Tm(x,  0)  —  9*)PmT{  •,  0)  dx. 


VmTmVPmfdx  <  \\vm  ||  L2  \\Tm  ||  L3  || VPmf  ||  L6  <  C\\vm  ||  \\Tm  ||: 1/2  \\Tm  ||  %  ||  PmT  ||  H2  . 


It  is  observed  that 


dPmT  PmT 

— - —  =  —a - ,  at  2  =  (J , 

OZ  H2 

PmT  =  0  ,  on  OM0  x  [-H0, 0] , 


dPmT 


=  0 ,  at  2:  =  —  H0 
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since  PmT  =  J2JLi  Pj^j  and  each  satisfies  the  boundary  condition  in  (3.6).  Consequently,  an  application 
of  the  elliptic  regularity  for  PmT  gives 


(3.21) 


\Pmf\\H2  <C\\A(PmT)\\. 


In  more  detail,  we  have 

m  m  m 

A(Pmf )  =  a(JT  ft*,)  =  J2  PiA*j  -  E  ’  by  (3-6)  ’ 


(3.22) 


3= 1  3  = 1  3= 1 

oo  oo  oo 

AT  =  a(Yj  =  E  =  E  Vl*.  >  by  (3.6) , 


3  =  1 


3= 1 


J=1 


which  along  with  the  orthogonality  of  {<I)J}?>1  in  L2(A4)  leads  to 


(3.23) 


\\A(PmT) f  =  E  A?42 Ill'll 2  <  E  Al^2 Ill'll 2  =  11^ II5 

i=i  j=i 


The  combination  of  (3.21)  and  (3.22)  results  in 

(3.24)  \\Pmf\\H-  <  C\\AT\\  <  C\\f\\H2  . 


The  substitution  of  (3.24)  into  (3.19)  leads  to 


(3.25) 


VmTm  VPrnT  dx 


' M 


<  C||um||||Tm||b'12||Tm||1/2||T||H2  <  C||Tm||^||Tm||1/2||T||ff2  . 


3/2ll T  II1/2! 


Similarly,  we  have  the  following  estimates 

Wrn  J-  m  @ zPmP  dx 


(3.26) 


IM 


<  C'||Tm||3ff/12||Tm||1/2||f||H2, 


(3.27) 


(vTm,VPmr)  +  (dzTm,dzPmf 


<  C\\Tm\\H4f\\H2 , 


(3.28) 


rv 

—  (Tm(x,0)-9*)PmT(;0)dx 

«2 


< 


C\\Tm\\Hi\\f\\H2 , 


By  the  combination  of  (3.16),  (3.17),  (3.25)-(3.28),  we  arrive  at 

(3.29)  dtTm  €  a  bounded  set  of  L4//3  ^0,  t\ ;  , 

for  any  t\  >  0  and  independent  of  m. 
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The  estimates  (3.16),  (3.17),  (3.29)  imply  the  existence  of  T  £  L°°(0,  t\ ;  L2)flL2(0,ti ;  H1)  and  u,v  ,w  G 
L2(0,ti ;  L'2)  and  a  subsequence  {Tm>,  vmi,  wm>}  such  that 

Tm>  T  weakly  in  L2(0,  t\ ;  H1) , 

(3.30)  Tmf  T  weak-star  in  L°°(0,  t\ ;  A2) , 

vm> ,  Wm'  — v  ,  w  weakly  in  A2( 0,  t\ ;  L 2) , 

With  the  use  of  (3.29),  (3.30)  and  Aubin’s  compactness  theorem,  we  also  have 

(3.31)  Tm /  — >  T  strongly  in  L2(0,  t\  ;  L2) . 


Then  it  is  standard  to  pass  to  the  limit  in  (3.9)-(3.11)  and  prove  that  the  limit  function  (T,v,w)  is 
indeed  a  weak  solution  as  defined  in  (3.2).  The  details  are  omitted  for  brevity.  The  proof  for  the  first  part 
of  Theorem  3.1  is  completed. 

Furthermore,  we  multiply  (3.15)  by  4>(t),  where  (f>  £  <j>(t)  >  0,  and  integrate  in  time: 


(3.32)  J  1  (||Tm(-,  <)||2  +  2  J  (ki|| VTm(-,  s)||2  +  K2\\dzTm(;  s)||2)  ds)  </>(t)  dt  <  j*  (||T0||2  +  C*t)</>(t)  dt . 


Using  the  weak  convergence  (3.30)  we  pass  the  lower  limit  in  this  inequality  and  obtain 

(3.33)  J*1  (||T(-,  t)||2  +  2  j*  (^11  VT(-,  s)||2  +  k2|| dzT(-,  s)||2)  da)  dt  <  £  (||T0||2  +  C*t)cj>(t)  dt , 

for  all  <j)  £  V^(0,ti)J ,  <j)>0.  This  amounts  to  saying  that  the  energy  inequality  (3.4)  is  satisfied  for  almost 
every  t  £  [0,  ti]. 

For  the  second  part,  we  note  that,  by  a  direct  calculation  using  the  representation  formulas  for  the 
horizontal  velocity  field  in  (3.2b), 

(3.34)  /  v(x,y,z)dz  =  ve(x,y) ,  V(x,  y)  £  Mo . 

J-H0 


This  leads  to  the  first  identity  of  (3.5),  due  to  the  free  divergence  of  ve  given  by  formulas  (2.10)  and  (2.13). 
Moreover,  the  second  identity  of  (3.5)  is  valid  since  ve  satisfies  the  specified  boundary  condition  on  the 
lateral  boundary.  The  third  identity  of  (3.5)  is  also  found  by  direct  verification  using  formula  (3.2c).  This 
completes  the  proof  of  Theorem  3.1.  | 
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